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Abstract 

We introduce the notion of a network's conduciveness, a probabilis- 
tically interpretable measure of how the network's structure allows it to 
be conducive to roaming agents, in certain conditions, from one portion 
of the network to another. We exemplify its use through an application 
to the two problems in combinatorial optimization that, given an undi- 
rected graph, ask that its so-called chromatic and independence numbers 
be found. Though NP-hard, when solved on sequences of expanding ran- 
dom graphs there appear marked transitions at which optimal solutions 
can be obtained substantially more easily than right before them. We 
demonstrate that these phenomena can be understood by resorting to the 
network that represents the solution space of the problems for each graph 
and examining its conduciveness between the non-optimal solutions and 
the optimal ones. At the said transitions, this network becomes strikingly 
more conducive in the direction of the optimal solutions than it was just 
before them, while at the same time becoming less conducive in the oppo- 
site direction. We believe that, besides becoming useful also in other areas 
in which network theory has a role to play, network conduciveness may be- 
come instrumental in helping clarify further issues related to NP-hardness 
that remain poorly understood. 

PACS numbers: 89.20.-a, 89.70.-a, 89.75.-k. 

1 Introduction 

The past decade has seen an impressive growth in the science of complex net- 
works, understood as the branch of scientific inquiry which, by merging well 
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established notions and techniques from the theory of graphs and from sta- 
tistical physics, addresses the interplay of structure and function in the large, 
essentially unstructured networks that occur in a wide variety of domains. The 
latter have encompassed several instances in many biological, social, and tech- 
nological fields, and have yielded an equally variegated array of results that the 
reader can now refer to in books and paper collections such as [51 HU [B] . 

One common methodological thread in all these studies has been the def- 
inition of a graph to represent the interactions among certain entities in the 
domain of interest, followed by the analysis of mathematical descriptors of some 
of the graph's properties as averages over a number of graphs generated ac- 
cording to some random-graph model thought to represent the phenomenon 
under consideration. Thus have emerged important finds regarding some net- 
works' characterization as small-world structures, or as scale-free structures, as 
well as powerful structural indicators of a graph's nature, such as its clustering 
coefficient and various centrality-related quantities. 

Here we introduce another indicator of a graph's properties, called its con- 
duciveness. Given a directed graph D of node set N, the conduciveness of D is 
defined with respect to two subsets A and B of N. Let di denote the out-degree 
of node i in D (the node's number of outgoing edges) and define I's B-bound 
out-degree, denoted by df , to be its number of outgoing edges whose heads 
are members of B (i.e., edges that lead from i to some member of B). The 
conduciveness of D from A to i? is denoted by Cond-A. b{D) and given by 

Cond^,s(I?) = ^'^-^5 . (1) 

Clearly, < Cond^,B(i:>) < 1. 

This definition of a directed graph's conduciveness can be easily interpreted 
in the context of hypothetical agents inhabiting the graph at its nodes but free 
to roam to other nodes by taking steps that follow edges along their directions. 
Specifically, Condyles (D) is the probability that, conditioned on there being one 
agent at each and every node in A, one further random step out of all those 
that they can take leads to a node in B. A graph for which this probability is 
higher than for another is regarded as more conducive to set B from set A than 
that other graph. 

Our initial motivation for the introduction of this definition has been its 
potential application to explain some phenomena related to the complexity of 
solving certain problems of combinatorial optimization. Normally such a prob- 
lem is defined on the set ft of the feasible solutions to the problem, using a real 
function / defined over Q according to which an optimal member of fl is to be 
found (one for which / is minimum or maximum over all of 51, depending on 
the problem). Many such problems are NP-hard, meaning that finding optimal 
solutions to them is at least as hard as solving any of the decision problems 
that constitute the class NP (those whose solutions, should they somehow be 
provided at no cost and turn out to be affirmative, could be checked to be cor- 
rect in polynomial time |14]V That no polynomial-time deterministic algorithm 
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has ever been found to solve an NP-hard problem is normally taken as a sign of 
computational intractability as problem instances grow large. 

This, however, is to be taken with caution. The class NP can be viewed as 
a complex hierarchy of subclasses |17j . which may ultimately help account for 
what is observed in practice: some NP-hard problems are solvable much more 
efficiently than others; more strikingly, two similarly sized instances of the same 
NP-hard problem may require considerably different amounts of computational 
effort to be solved. A way to illustrate this that is useful in the context of this 
paper is based on the following. Let there be n nodes and, for M — n(n — l)/2, 
consider a sequence Q — {G\,Gi: ■ ■ ■ , Gm) of undirected graphs on these nodes. 
For m = 1, Gm has one single edge joining two randomly chosen nodes and 
n — 2 isolated nodes. For 1 < m < M, Gm is obtained from Gm-i by placing a 
further edge between two randomly chosen nodes that are not already joined by 
an edge; so Gm has m edges and, for relatively small m, may also have isolated 
nodes. 

The crucial observation is that, as first documented in [521 H] in the wake 
of what was done earlier for some NP-hard decision problems [TT] [1^1 [201 [TSl 
[12], there exist NP-hard optimization problems for which practically every 
attempted algorithm, deterministic or otherwise, undergoes sharp transitions 
when applied to the graphs in Q for increasing values of to. These transitions 
refer to how long it takes the algorithm to reach an optimal solution and happen 
at well-defined values of m given Q. The same initial reporters of these phe- 
nomena also offered tentative explanations related to the nature and structure 
of the corresponding sets (one for each of the M graphs in ^), but those have 
lacked full consistency owing to normalization difhculties as the sizes of those 
sets grow along with to [4]. It has also been a difficulty that the values of to 
at which the transitions occur tend to be different if Q is changed, so the afore- 
mentioned analyses have only addressed single graph sequences and therefore 
lack statistical significance as well. 

We have found that the notion of a directed graph's conduciveness, as intro- 
duced above, has an important role to play in elucidating the nature of these 
optimization transitions. The fundamental idea is, for each to, first to identify 
an appropriate descriptor of the feasible solutions to the problem that is being 
posed on Gm- This will give us the VL set for that particular m, henceforth de- 
noted by f2m. Then we identify some primitive operation on the members of VLm 
that may be used to transform one of them into another. Every two members of 
rim that are thus related constitute an ordered pair; collectively, all such pairs 
constitute a set that we denote by Em- The directed graph whose conducive- 
ness we study, denoted by Dm, has node set ^Im and edge set Em- This graph 
embodies all primitive steps that an optimum-seeking algorithm may take to 
solve the problem on Gm- For m = 1,2,..., M, we study the conduciveness of 
Dm from the nodes in Q.m that do not represent optimal solutions to those that 
do. 

By its very nature as a probability, this conduciveness of Dm has none of 
the normalization problems alluded to above. As we will see, it also allows for 
some multiplicity of events to be investigated for statistical significance, though 
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only to a limited extent. This is because in general the edges of can only be 
found through the explicit enumeration and testing of several pairs of members 
of firm which in most cases is a very large set even for very small values of n. 
We then see that there exist severe time constraints on the generation of the Dm 
graphs for multiple instances of the sequence Q, and consequently constraints 
on the largest value of n that can realistically be used. We note, moreover, 
that seldom can Dm be fully stored, which limits the properties that can be 
analyzed. 

We target the same two optimization problems as [3], namely the prob- 
lem of coloring the nodes of an undirected graph optimally and that of finding 
one of its maximum independent sets. Aside from the fact that they are both 
paradigmatic NP-hard optimization problems, our choice of them has also been 
influenced by the remarkable fact that, for each value of to, it is possible to define 
a single rim set for both problems, thus allowing the study of their optimization 
transitions to be conducted in a peculiarly interrelated fashion. 

2 Graph coloring and independent sets 

The chromatic number of an undirected graph G on n nodes, denoted by x(G), 
is an integer between 1 and n indicating the smallest number of distinct colors 
(labels) that can be used to tag the nodes of G in such a way that every node 
gets exactly one color and no two nodes connected by an edge get the same color. 
The independence (or stability) number of G, denoted by a(G), is likewise an 
integer between 1 and n and indicates the size of a largest independent subset 
of G's node set, that is, a largest subset of nodes containing no two nodes 
connected by an edge [7]. Finding either number is an NP-hard problem |14j . 

The two problems can be reformulated in such a way that their sets of 
feasible solutions are in fact the same set. To see this, first let an orientation of 
G be an assignment of directions to G's edges, that is, one of the ways in which 
G can be turned into a directed graph. An orientation is acyclic if it contains 
no directed cycles (i.e., it is never possible to return to a node after moving 
away from it along the directions of the edges). Every acyclic orientation of G 
yields a number of colors to tag the nodes of G legitimately, and likewise an 
independent set of G. Conversely, every legitimate assignment of colors to the 
nodes of G yields an acyclic orientation of G, and so does every independent 
set of G. The proofs that back up these statements are not simple [13], but 
accepting them clearly implies that both finding x(C) ^'^d finding a(G) can 
be formulated based on sharing the 51 set defined as the set of all the acyclic 
orientations of G. 

The precise relationships implied by the proofs in il^] are the following. Let 
uj be an acyclic orientation of G. Let Depth(w) denote the number of nodes on 
a longest directed path in G according to uj, henceforth referred to simply as 
the depth of uj. Then 

X(G) =minDepth(c^), (2) 
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that is, x{G) is the depth of the shahowest member of $7. Now let Width(w) 
denote the least number of node-disjoint directed paths into which G can be 
decomposed given w, henceforth referred to simply as the width of uj. Then 

a(G) =maxWidth(a;), (3) 

meaning that a{G) is the width of the widest member of VI. 

It also emerges from those same proofs (but see [2, 3 for explicit accounts 
of the corresponding algorithms) that, given both Depth(w) and Width(a;) 
can be computed in polynomial time. So, by Eqs. ^ and the NP-hardness 
of the two problems in question is to be attributed to the inherent difficulty 
of searching inside for an optimal uj in each case. Following the general 
outline provided in the previous section, we continue our analysis by defining 
the directed graph D of node set D, whose edge set, is to be set up to 
reflect some primitive relationship among the members of 51 that can be used 
to transform each one into some other. 

There are certainly several ways in which an acyclic orientation, say cj, can 
be turned into another, say uj' . One possibility that has become popular in 
several task scheduling applications is to turn one or more of the sinks of w 
(nodes with no outgoing edges) into sources (nodes with no incoming edges) 
and then let the resulting orientation be oj' (clearly acyclic, given the acyclicity 
of w). We eschew this choice for two reasons. The first one is that it entails 
several direction reversals for one single transformation, which then seems hard 
to qualify as primitive. The second reason is that, under such sink-to-source 
transformations, the resulting D is almost always a fragmented graph (i.e., there 
exist pairs of nodes that are unreachable from each other even if edge directions 
are ignored) [T] . Since our interest is in the conduciveness of D with respect to 
certain subsets of J7, it seems that starting out with a fragmented D is bound 
to produce results somewhat devoid of meaning. 

Our definition of the edge set i? of D is then the following. Given cj G fl, 
an edge exists directed from u to some cj' £ if and only if uj' results from 
reversing the direction of exactly one of the edges of G as oriented by cj. Of 
course, if (w,a;') £ E holds, then so does (w',aj) G E. Moreover, it now holds 
that D is strongly connected (that is, a directed path exists from any node to 
any other) 

Handling D computationally, though, is a difficult matter owing to both its 
number of nodes and the explicit way in which its edges must be enumerated. 
The number of nodes, which is the number of distinct acyclic orientations of 
G, is given by a surprising application of the so-called chromatic polynomial of 
G [23| and, for n fixed, grows rapidly from the two orientations allowed by the 

^If not, then there have to exist iO,ui' S f2 with the following property. If S is the set of 
edges in G on whose directions uj and uj' disagree, then the individual direction reversal of any 
edge e a S creates a directed cycle C in the resulting orientation. But since cj' is acyclic, C 
must also contain another of the edges of S, say e', and this edge's direction must oppose that 
of e on C. Notice, however, that C comprises at least three edges, so reversing the direction of 
e' alone would create no directed cycle. This contradicts the existence of the uj, u)' pair with 
the assumed property. 
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case of one single edge to the n\ orientations that a graph with all possible M 
edges on n nodes admits. As for discovering the edges of D that outgo from a 
particular orientation w, there is in general no alternative but to try and reverse 
the directions of all edges of G, one by one with respect to what w stipulates, 
checking for each one whether the resulting orientation is itself a member of fi. 

Given G, we enumerate the members of fl by the algorithm given in but 
store each one only while recording some of its properties for later use. For 
each oj that is output by the algorithm, we calculate Depth(w), Width(a;), its 
out-degree di^ in D, and its S-bound out-degree d^- Here B depends on which 
problem is being addressed. If it is the coloring problem, then B is the subset 
of comprising orientations whose depths are all equal to x(G). If it is the 
independent-set problem, then B is the subset of fl whose orientations all have 
width a{G). For simplicity, whenever the context allows we refer to di^ as a full 
degree and to dj^ as an optimum-bound degree. Note that each full degree is an 
integer between 1 and the number of edges of G. An optimum-bound degree, 
in turn, is an integer between and again the number of edges of G. 

Scarce though they may be, these recorded properties of D allow for some 
useful statistics to be computed, in addition to allowing for the direct calcula- 
tion of some useful conduciveness figures for D as per Eq. Using 6x.y to 
denote Kronecker's delta function of the integers x and y, and \X\ to denote the 
cardinality of set X, these statistics are: 

• The distribution of full degrees in Z?, given by 

P{k) = ^^^gp^ (4) 

for every possible full degree k. 

• The distribution of depths in D, given by 



Q{u) 




,U 



(5) 



for every possible depth u. 



• The distribution of widths in D, given by 



R{v) 




,V 



(6) 



for every possible width v. 



• The joint distribution of full degrees and depths in D, given by 



S{k, u) 




,u 



(7) 



for every possible full degree k and depth u. 
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• The joint distribution of full degrees and widths in D, given by 

^ - w 

for every possible full degree k and width v. 

Additional statistics are P^{k), Q^{u), R^{v), S^{k,u), and T^{k,v), defined 
analogously to the above but for optimum-bound degrees (that is, substituting 
for in the corresponding definitions). Note also that, whenever the G 
in question is one of the Gm graphs of the sequence Q introduced previously, 
we alter the notation of these statistics by adopting the subscript m for them 
as well (consistently with the graph Dm of node set and edge set Em, all 
introduced earlier but now with the specific meanings given in this section for 
D, f], and E, respectively). 



3 Optimization transitions in random graphs 

Let us then look at one single sequence Q — (Gi, G2, . . . , G^/) for n ~ 12 
(whence M = 66) and observe how different algorithms to find x(G'i), xiGi), ■ ■ ■ , 
xiCrAi) and a(Gi), a(G2), . . . , a(GM) perform. Results on finding the chromatic 
numbers are given in Fig. [1] those on finding the independence numbers are in 
Fig.H 

Figure [1] contains performance data on three algorithms. First is a simple 
random walker, which for each m starts at a randomly chosen acyclic orientation 
in Qm and then at each step traverses one of the edges that outgo from its current 
acyclic orientation in the set Em, thus reaching another acyclic orientation. 
Performance data are then given for a genetic algorithm operating on ilm, and 
then for a deterministic algorithm whose operation is not based on Dm at all. 
The random walker and the genetic algorithm stop upon hitting the first acyclic 
orientation uj for which Dcpth(w) = x(G,„), with the provision that x(Gm) 
is known beforehand from running the third algorithm first. The data on all 
three algorithms are shown in the three parts of Fig. [1] against a backdrop of 
vertical lines, each marking the number of edges at which a transition in the 
chromatic number occurs: for 1 < to < M, if x{Gm+i) = x(G'm) + Ij then a 
vertical line is drawn at the abscissa to + 0.5. The chromatic numbers of the 
graphs in Q necessarily increase by 1 at each transition, from x{Gi) = 2 through 
x{Gm) = n, therefore there are n — 2 vertical lines all told. 

A similar arrangement holds for Fig. [51 whose setting differs from that of 
the previous one in that now both the random walker and the genetic algorithm 
stop upon finding uj 6 flm such that Width(w) = a{Gm), once again given 
that a{Gm) is known a priori from running the deterministic algorithm first. 
There is also an important difference regarding the transitions in the graphs' 
independence numbers, which now necessarily decrease by 1 at each transition, 
from a(Gi) = n — 1 through a(GM) = 1, thus totaling n — 2 transitions as 
well. So, in Fig. [2l the vertical lines marking the transitions are drawn at the 
abscissae m + 0.5 such that a(Gm+i) = a{Gm) — 1. 
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Parts (a) and (b) of both Figs. [T] and [5] thus refer to methods which work 
on the sets fim of acychc orientations of the graphs Gm , either making exphcit 
use of the structure of each (which the random walker does) or ahowing for 
longer jumps as the orientations undergo the crossover and mutation operations 
prescribed by the genetic algorithm. The data displayed in the corresponding 
four panels often have in common the property that the occurrence of a transi- 
tion, say from m' to to' + 1 edges, causes the algorithm in question to perform 
significantly better on Gm'+i than on Gm' , and then increasingly poorly through 
the following values of m until the next transition, if any, is reached. This dif- 
ference in performance is sometimes quite marked, involving improvements by 
at least one order of magnitude. 

What is perhaps more curious is that the same behavior is also present in 
Fig. [ijc) , which refers to a deterministic algorithm to find chromatic numbers 
that does not rely on the Dm graphs (in fact, this algorithm's underlying strat- 
egy makes no reference at all to the acyclic orientations of the graph whose 
chromatic number it is seeking). Informally, then, this seems to indicate that 
the characteristic performance jumps at the transitions are inherent to the op- 
timization problem itself (and only marginally, if at all, dependent upon how 
its feasible solutions are represented). It also seems to confer to the Dm graphs 
some of the primitive representational character we sought in the beginning. 
However, Fig. HJc), which also refers to a deterministic algorithm that does 
not operate on acyclic orientations, only now to find the graph's independence 
number, shows none of the effects on performance at the transitions that the 
random-walk and genetic-algorithm approaches exhibit. The reason for this is 
that, despite being just as nominally NP-hard as the problem of finding chro- 
matic numbers, finding independence numbers is easier in practice than that 
other problem. What this means is that, in order for the transitions' effects on 
performance to show, substantially higher values of n are needed (cf. Fig. 9 in 
i)- 

We then proceed on the premise that the sequence of Dm graphs for to = 
1, 2, . . . , M contains information enough to explain the performance jumps at the 
transitions, even though quantitatively we can only resort to what can be derived 
from each graph's nodes' depths, widths, and out-degrees. One initial indication 
that this makes sense comes from investigating the mutual information of pairs of 
random variables associated with each Dm- Given two discrete random variables 
and the joint distribution of their values, their mutual information is a measure 
of how much fixing the value of one of them reduces the uncertainty on the value 
of the other |18| . Just like Shannon's entropy, mutual information is expressed 
in (information-theoretic) bits. 

In the context of finding x(G'm), two discrete random variables of interest 
are those that give the out-degree and the depth of a randomly chosen node of 
Dm- Their joint distribution is Sm{k,u) in the case of full degrees, Sm{k,u) 
in the case of optimum-bound degrees, both introduced earlier. Their mutual 
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information is, respectively for each case, given by 



k=l u=2 

and 



ZZSS.ik.ui"'l,^i^y (10) 



/c=0m=2 



If the problem is to find a{Gm), then the two random variables give the node's 
out-degree and its width. Their joint distributions are T,„(fc,?;) and T^{k,v), 
once again depending on whether full or optimum-bound degrees are referred 
to. The respective measures of mutual information are 

m n-1 T (h \ 



and 

jS,-f:t^S,ik..,)V.,.j§^y (12) 



A:=0 11=1 

For the same sequence Q we have considered so far in this section, and follow- 
ing the same conventions as Figs. [T] and [5] with regard to marking with vertical 
lines the values of m at which x{Gm) or a{Gm) changes along the sequence, we 
show in Fig. [3] the progress of these four mutual information functions, in part 
(a) for graph coloring, in part (b) for independent sets. Note, in all cases, that 
although consistently less than one bit, all four functions are nearly always pos- 
itive, thus providing evidence that, in most Dm instances, out-degrees of either 
kind are independent from neither depths nor widths. However, the functions 
do not seem to behave consistently at the transitions and for this reason offer 
no direct explanation for what happens there. 

It is important to recall that all of Figs. [T]-[3] refer to the one single sequence 
Q. They have been offered as illustrations of what is typical, and averaging over 
multiple graph sequences, which requires that we address the fact that a given 
transition may occur at different m values in different sequences, might blur the 
reader's understanding of what the phenomenon is and how mutual information 
suggests that it has to do with the graphs. We now turn to the role played 
by graph conduciveness and, after one more single-sequence illustration, do some 
averaging as properly as possible. 



4 Computational results on conduciveness 

Given a Dm graph, let denote the subset of Vim whose members are those 
of depth x{Gm)- Analogously, let fi^ denote the subset of ^Im whose members 
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are those of width a{Gm)- We study four kinds of conduciveness of D,n, given 
as follows with reference to Eq. (IT|)o 

C^'- = Condn„\a-,n-(^™); (13) 

C^°^' = Condax,^^\^x(i?„0; (14) 

C^'^ = Conda„\no,n=(i?„0; (15) 

C;^-™' = Conda=,a„An=(^m)- (16) 

Note that C^'™ is the conduciveness of Dm from all orientations that are non- 
optimal for coloring to those that are optimal, C^'°"* the conduciveness in the 
opposite direction. The situation with C^^'" and C^'°"* is totally analogous, 
now regarding optimality for independent sets. We refer to C^''" and C^''" as 
being inbound, to C^'°"* and C^'°"* as being outbound. Note also that, by 
Eq. ([T]) , and given the antiparallel nature of the edge set of , it holds that 

^XiOut ^cj£Sl„i\0^ ^ ^xA^^ (17) 

and 

<-^m — — -1 ■ \^°) 

These, however, imply no obvious relationship between the inbound conducive- 
ness of Dm and its outbound conduciveness for any of the two problems. 

An illustration of the kind of relationship that does hold is given in Fig. 21 
which results from our last use of the same single sequence Q as in the previous 
section. Strikingly, as the sequence unfolds with increasing m and the transi- 
tions in x{Gm) [part (a) of the figure] or a{Gm) [part (b) of the figure] occur, 
the inbound conduciveness of Dm undergoes sudden jumps upwards precisely at 
the transitions while its outbound conduciveness undergoes downward jumps. 
The inbound-conduciveness jumps can be seen to encompass at least one order 
of magnitude in many cases. Between one transition and the next, the inbound 
conduciveness deteriorates progressively while the outbound conduciveness im- 
proves. This is then the key to interpreting the phenomena illustrated in Figs.[T] 
and[2l at the transitions. Dm becomes markedly more conducive in the direction 
of the optimal orientations and less conducive in the opposite direction; right 
past a transition through right before the next one happens. Dm tends to be- 
come progressively less conducive in the direction of the optima, more conducive 
in the direction that leads away from them. 

Next we study these conduciveness variations as averages over the graph se- 
quences Gi,G2, ■ ■ ■ ^ Gi5, each comprising graphs on n = 12 nodes and generated 
independently. As noted earlier, even though both the chromatic number and 
the independence number undergo n ~ 2 transitions each in each sequence, the 
tth transition, for some iG{l,2,...,n — 2}, may happen at different values of 



\ denotes set difference. 
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m for the different sequences. Some alignment of the transitions is then needed 
for the averages of interest to be computed; we proceed as follows. If for a 
given sequence the tth chromatic-number transition occurs for m = m! , then we 
calculate the change ratios iC^, - C^)/C^ and {C^?l\ - 
We do likewise for each independence-number transition. If two subsequent 
chromatic- number transitions occur at m — m' and m — m" > m' , then we also 
calculate the change ratios and {C^J'^' - C^J^) / C^J^ 

again proceeding likewise for the interval between every pair of subsequent 
independence-number transitions. The latter formulae can also be used to cal- 
culate change ratios for the interval that precedes the first transition (letting 
m' = and m" be the value of m at which the first transition occurs) and the 
interval that succeeds the last transition (letting m' be value of m at which the 
last transition occurs and m" = M). Once all change ratios have been calcu- 
lated, they can be averaged over the 15 sequences for each transition (whichever 
the value of m is at which it happens to occur in each sequence) and each 
interval. 

These average change ratios are given in Figs.[5]and|ni respectively for graph 
coloring and independent sets. All data in these two figures are presented, as 
in all previous cases, against a backdrop of vertical lines. These, however, are 
now equally spaced and refer to the transition numbers, from 1 through n — 2, 
regardless of the m values at which the transitions themselves are observed in 
each particular sequence for each problem. Each panel in each figure contains 
two plots, one with points whose abscissae coincide with those of the vertical 
lines (this refers to change ratios at the transitions) and one with points whose 
abscissae stand either halfway between those of two consecutive vertical lines or 
to the left (right) of the leftmost (rightmmost) vertical line's abscissa [this refers 
to change ratios along the intervals between consecutive transitions or before 
(after) the first (last) transition]. 

Figs.[5]Ja) andlUa), which refer to the progress of the inbound conduciveness 
as the transitions elapse, reveal that at most transitions the upward jumps 
represent significant fractions of the pre-transition conduciveness values, which 
often increase manyfold (by a factor of a few tens). On the other hand, the 
accumulated deterioration in conduciveness that is observed between transitions 
and in the outermost intervals is in most cases given by a fraction that varies 
widely depending on the transition, ranging practically from nearly no loss of 
the initial conduciveness inside the interval to nearly total loss. 

Figs. El^b) and[6l^b), in turn, refer to how the outbound conduciveness val- 
ues evolve along with the transitions and show that, right at the transitions, 
conduciveness is lost with respect to the pre-transition values by fractions that 
amount to losing from about 5-10% of it (depending on the problem) to all 
of it. As we look at the accumulated improvement in conduciveness between 
transitions and in the outermost intervals, we see that a wide range of possibil- 
ities is again present, allowing at one extreme for practically no improvement 
and, at the other extreme, for an improvement by about 70-90% of the initial 
conduciveness inside the interval (depending on the problem). 
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5 Summary and outlook 



The notion of network conduciveness we have introduced is a simple degree- 
based indicator that can be interpreted as a probabihty with respect to a par- 
ticular agcnt-rchitcd dynamics. We believe that, either as defined or as some 
variant thereof, it may find applications in network studies having to do with 
the dynamics of populations in networks. Our own application in this paper 
has been to the field of combinatorial optimization, and then the network in 
question is representative of the feasible solutions to a particular instance of an 
optimization problem and of how one may move from one sohition to another 
through as simple a local transformation as possible. We tackled the NP-hard 
problems of finding an undirected graph's chromatic and independence numbers 
and demonstrated how network conduciveness, when applied to problem repre- 
sentations in the domain of the graph's acyclic orientations, is capable of helping 
explain the well-known performance transitions that occur along sequences of 
random graphs for both problems. 

As it happens, though, the networks whose conduciveness we have considered 
grow very rapidly with the graph's numbers of nodes and edges, and become 
themselves very nearly intractable already for small instances of the problems 
we addressed. We were then limited in our computational experiments to using 
graphs on 12 nodes exclusively and to averaging results on 15 sequences of ran- 
dom graphs. For the sake of the record, with current technology all experiments 
required nearly two months on twenty processors. So, as much as we think that 
there is great potential usefulness to the notion of a network's conduciveness, 
further progress with the particular application we chose requires considerable 
further effort so that larger graphs and better statistical significance can be 
aimed at. On the other hand, we regard the first steps we have taken as very 
significant: to the best of our knowledge, no other study has addressed the intri- 
cacies of NP-hard optimization problems from the perspective of network theory 
applied to the structure that underlies the problems' sets of feasible solutions. 
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Figure 1: Performance of three algorithms to find the chromatic numbers of the 
graphs in ^: a random walker (a), a genetic algorithm (b), and Trick's imple- 
mentation [24] of the algorithm in [9] (c) . The data shown for (a) are averages 
over 10 000 independent runs. The genetic algorithm is the one described in the 
caption to Fig. 15 in [T, parameter values included, and the data shown for (b) 
are averages over 100 independent runs. 
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Figure 2: Performance of three algorithms to find the independence numbers of 
the graphs in ^: a random walker (a), a genetic algorithm (b), and Johnson's 
implementation [16] of the algorithm in jTO] (c). The data shown for (a) are av- 
erages over 10 000 independent runs. The genetic algorithm is the one described 
in the caption to Fig. 15 in [5, parameter values included, and the data shown 
for (b) are averages over 100 independent runs. 
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Figure 4: C^"' and (7^'°"* (a); C^''" and C^'°"* (b). All data refer to the 
sequence Q. 
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Transition 

Figure 5: Conduciveness change ratios at the graph-coloring transitions and 
related intervals. Data are given as averages over the set of sequences 
Gi,Q2t ■ ■ 1 Gi5 for both the inbound conduciveness (a) and the outbound con- 
duciveness (b). 
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Figure 6: Conduciveness change ratios at the independent-set transitions 
and related intervals. Data are given as averages over the set of sequences 
Gi,Q2t ■ ■ 1 Gi5 for both the inbound conduciveness (a) and the outbound con- 
duciveness (b). 
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